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ABSTRACT 

Gamma-ray bursts are believed to originate from internal shocks which arise 
in an irregular relativistic wind. The process has been thought to be inefficient, 
converting only a few percent of the kinetic energy into gamma-rays. We define 
ultra efficient internal shocks as those in which the fraction of emitted energy 
is larger than the fraction of energy given to the radiating electrons at each 
collision. We show that such a scenario is possible and even plausible. In our 
model, colliding shells which do not emit all their internal energy are reflected 
from each other and it causes subsequent collisions, allowing more energy to 
be emitted. As an example, we obtain about 60% overall efficiency even if the 
fraction of energy that goes to electrons is e e = 0.1 provided that the shells' 
Lorentz factor varies between 10 and 10 4 . The numerical temporal profile 
reflects well the activity of the source which ejects the shells, though numerous 
collisions take place in this model. 



Subject headings: gamma rays: bursts; shock waves; relativity 



Introduction 



A widely accepted mechanism for producing a cosmological gamma-ray burst (GRB) 
is the deceleration of relativistically expanding shells. The kinetic energy of the shells 
is converted into internal energy by relativistic shocks. These shocks can be due to 
collision with the ambient medium (external shocks) or shocks inside the shell itself due to 
nonuniform velocity (internal shocks). Electrons are heated by the shocks, and the internal 
energy is radiated via synchrotron and IC emission, with broken power law spectra (e.g. 
Sari, Piran and Narayan 1998). 

Most bursts have a highly variable temporal profile with variability time scale 
significantly shorter than the overall duration. In the external shocks scenario this 
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variability is due to irregularity in the surrounding material, but the efficiency is extremely 
low (Fenimore, Madras & Nayakshin 1996; Sari & Piran 1997). Thus, GRBs are believed 
to be produced by the other alternative: internal shocks. The inner engine itself should be 
variable in this scenario, the observed temporal profile follows very closely the operation 
of the source (Kobayashi, Piran & Sari 1997, KPS97 hereafter). Using a simple model, 
we have estimated that the hydrodynamic efficiency of this process (transforming kinetic 
energy to internal energy) is about 10% . 

Kumar (1999) argued that the conversion efficiency from the bulk motion to gamma-ray 
is only 1%. His argument is based on three points: (I) The hydrodynamic efficiency is, as 
mentioned above, typically 10%. (II) It is only the electrons that are radiating. Even in 
equipartition among protons, magnetic field and electrons, the latter electrons only have a 
third of the internal energy. (Ill) The amount of the radiated energy within the gamma-ray 
band is about a third of the total radiated one. Combining these three factors give the low 
efficiency of ^ x | x | = 1%. 

Such a low efficiency results in sever energy demands on the source. Moreover, it has 
problems explaining the energy ratio between the GRB and its afterglow. According to the 
internal-external shock model, the remaining kinetic energy, which was not converted to 
radiation by internal shocks, is radiated during the afterglow stage. External shock does 
not suffer from problem (I), and the energy released in the afterglow should be considerably 
higher than that in the GRBs. However, it seems to be that the energy during the afterglow 
is only a tenth of that during the GRB, rather than ten times larger (Frontera et. al. 
2000, Kumar & Piran 2000, Freedman & Waxman 2000). Even though the observational 
constraints are not very good since most of the energy is released at very early radiative 
stages where the afterglow was not observed, a factor of ten more energy in the afterglow 
seems to be excluded. 

A possible solution to these problems is to assume large angular fluctuation in the 
shells (Kumar and Piran 1999). This model has clear predictions in the form of afterglow 
variability whose amplitude decay in time. It also predicts that the afterglow may be 
sometimes more energetic than the GRBs. 

In this paper, we suggest a simple alternative solution, which overcomes the problems 
suggested by Kumar. We show that if the distribution of the Lorentz factor is not uniform, 
but instead its logarithm is distributed uniformly, then the typical ratio of Lorentz factors 
between neighboring shells is considerably larger. Then, the hydrodynamic efficiency can 
be close to hundred percents, even for a reasonable spread of Lorentz factors (Similar 
calculation was recently done by Beloborodov 2000). The main point of this paper is the 
possibility of "ultra-efficient" internal shocks. We define ultra-efficient internal shocks as 
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a scenario in which the emitted fraction of kinetic energy is larger than e e , the fraction of 
internal energy that is going into electrons (and then radiated) at each collision. We will 
show that such a scenario is possible and even reasonable. 

2. "Ultra Efficient" Internal Shocks 

Internal shocks could occur within a variable relativistic wind produced by a highly 
variable source. We represent the irregular wind by a succession of relativistic shells with 
a random distribution of Lorentz factors in a similar manner as in KPS97. Beloborodov 
(2000) has shown that the internal shocks can convert most of the kinetic energy to internal 
energy if the fluctuation of the initial Lorentz factors A 2 = ((7 2 ) — (l) 2 )/ (l) 2 is large. 
Though it is maximally 1/3 if the initial Lorentz factors take random values between 7 mm 
and 7 mm , it is not limited if the distribution is uniform in logarithmic space between 
log 7mm and logj max . 

The efficiency can be estimated by an equation similar to equation 19 in KPS97. The 
most efficient case is that the masses of the shells are taken equal, and that the efficiency is 
given by a simple form: 

(e)~l-a 1/a loga/(a-l) (1) 

where a = 'jmax/lmin- The efficiency is plotted as a function of the fluctuation 
A 2 = — 1 + (a + l)loga/2(a — 1) in figure [I]. This analytic estimate fits the result of 
our numerical simulation, and has the asymptotic form (e) = A 2 /2 which Beloborodov 
estimated for a small fluctuation. Our expression above generalizes Beloborodov estimate 
and gives reasonable estimate of the efficiency also for large fluctuations. 

Though we have seen that a large hydrodynamic efficiency is possible if the fluctuation 
of the initial Lorentz factors is large, it is not reasonable that all the internal energy is 
emitted after each collision, since electrons do not have most of the internal energy. Defining 
e e as the fraction of energy given to electrons we expect e e < 1. Even at equipartition 
with protons e e = 1/2. Under these circumstances the total emitted energy in Beloborodov 
model is still limited by e e . 

However, if e e < 1 the merger produced by a collision is expected to stay hot after the 
emission. As a result, the merger will spread to transform the remaining internal energy 
back to the kinetic energy (Kumar 1999). A simplified description of this process is to 
assume that the two shells reflect with a smaller relative velocity. The difference of the 
kinetic energy before and after the collision is the emitted internal energy. The reflecting 
shell will collide into the other neighbor shell. Since in this way a large amount of collisions 
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are caused, the overall efficiency from the kinetic energy to the radiation could be larger 
than e e . This is the key ingredient of our model. Before going on with this model, we 
present some hydrodynamic simulations which shows that these simplified assumptions are 
reasonable. 



3. Hydrodynamic Simulation 

To estimate the conversion efficiency of the internal shocks process, it is important to 
understand how high velocity shells interact with slower ones and dissipate the the kinetic 
energy. We consider here a collision of two equal mass shells with very different Lorentz 
factors, 7r/7s ^> 1. The rapid and the slower shells are denoted by the subscriptions r and 
s respectively. We assume that the widths of the shells are comparable in the ISM rest 
frame. This is a reasonable assumption since in this frame, the width is given directly by 
the "inner engine". Even so, the width of the rapid shell l r is much larger when viewed 
from the rest frame of the rapid shell l r /l s ~ (lr/ls) 2 /2. Then, the dense "slower" shell 
with Lorentz factor j s ~ i r /^ls collides into the low density "rapid" shell in this frame, 
Ps/Pr ~ lr/ls- This is a planar analog of the evolution of a relativistic fireball (Sari & Piran 
1995, Kobayashi, Piran and Sari 1999). When it begins to interact with the surrounding 
material, two shocks are formed: a forward shock propagating into the rapid shell and a 
reverse shock propagating into the slower one. After the reverse shock crosses the slower 
shell, the profile of the shocked rapid shell material approaches to that of its fireball analog: 
the "blast wave" which sweeps and collects the surrounding material. Once the shock wave 
crosses the rapid shell, the hydrodynamical structure is as follows. There is a shocked rapid 
shell, the analog of the "blast wave" and the shocked slow shell, which cooled down and is 
the analog of the adiabatically cooling "fireball ejecta" . The structure of such a system in 
the fireball case was studied by Sari and Piran 1999 and Kobayashi and Sari 2000. 

Since 7 2 3> p s / p r , the slow shell is considerably decelerated by the relativistic reverse 
shock down to ~ (lr/ls) 3 ^ and heated to a relativistic temperature at the crossing time 
t s ~ {ir/lsY^ hi c. Then, it cools adiabatically and follows a planar version of the Blandford 
Mckee (1976) solution (Sari 2000), in which a given fluid element evolves with a bulk Lorentz 
factor of 7 oc /~ 3 / 2 . Therefore, the motion becomes Newtonian 7 ~ (lr / ls) 3 ^ 4 (tr / ts)~ 3 ^ 2 ~ 1 
when the forward shock crosses the rapid shell at t r = l r /c. On the other hand, since the 
forward shock itself evolves as 7 oc Z^ 1 / 2 , it slows down to 7 ~ (^ r /^ s ) l l 2 /2 at the crossing 
time t r . 

We have developed a relativistic code with an exact Riemann solver to solve relativistic 
hydrodynamics problems (Kobayashi, Piran & Sari 1999). Using this code, we numerically 
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study the collision of two equal mass slabs with 7 S = 10, 7 r = 10 3 and the same width in 
the ISM frame. The initial condition in the rapid shell comoving frame is 

7r = 1 7s ~ 50 

p r = 1 p s = 100 

i r ~ 5000 i s = i 

The mass density outside the slabs and the homogeneous pressure are p = 10~ 8 and 
p = 10 -10 . The mass density and the pressure are measured in the comoving frame of 
each fluid. Our adiabatic simulations represent a case in which e e is very small. In this 
simulation, we used the unit of c = 1. 

To compare the analytical estimates with the numerical simulation we define the 
effective Lorentz factor of each shell as (7) = / vn^dlj J mdl with the effective mass 
m = 7{p + (3 + /3 2 )p}. The numerical simulation then gives (%) ~ 5.3 and (%) ~ 1.6 at 
the crossing time, which are in agreement with the analytical estimates. The thin line in 
figure H shows the numerical density profile at this time. After the forward shock crosses 
the rapid shell, a rarefaction wave begins to propagate into the rapid shell to transform the 
internal energy to the kinetic, then the rapid is accelerated to 7 ~ / ~f r /2 / ~f s . The center of 
mass moves with a Lorentz factor ~ ( , y r / r y s ) /2 ~ 5 in the rapid shell comoving frame. At 
the end of the simulation t = 10 5 , about forty percent of the rapid shell material is slower 
than the center of mass, and goes with the slower shell (see fig §b). 

4. Two Shell Collision 

A collision of two shells is the elementary process in our model. A rapid shell catches up 
a slower one and the two merge to temporarily form a single one (denoted by the subscript 
m). Using conservation of energy and momentum, the Lorentz factor of the merged shell 
7 m and the internal energy E int produced by the collision are given by 

7m ~ \J (mr'Jr + m s ^ s )/{m r /^ r + m s /7 s ) , E int = m r {j r - j m ) + m s {j s - j m ). (2) 

After a fraction e e of the internal energy is emitted isotropically in the local frame 
of merged shell (center of mass frame), the shells will spread to transform the remaining 
internal energy back to kinetic energy. If the widths are the same in center of mass frame, 
each mass is conserved before and after the collision. However, as we have seen, some 
fraction of the rapid shell material goes together with the slower one after the collision in 
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general. We parameterize the mass splitting as m' r = (1 — 8)m r and m' s = m s + 5 m r . We 
later show that the total efficiency is not sensitive to 5. 

The Lorentz factor of the reflected shells in center of mass frame are given by 

f r = (M 2 + m' r 2 - m' s 2 )/2m' r M, T s = (M 2 + m' s 2 - m' r 2 )/2m' s M. (3) 
where M = im r ^ T + m s '~f s — e e Ei nt / c 2 )/^ m . The Lorentz factors in laboratory frame are 

r r = f r7m - ^/(f2-i)(72,-i), r s = f slm + ^/(f2- 1)^-1). (4) 

The shells are once compressed by shocks, but these spread when reflecting. For simplicity 
we assume that the width of the shells, k, is constant. 



5. Multiple Shell Collision 

We consider a wind consisting of iV shells. Each shell is characterized by the four 
variables: / -f i ,m i ,l i and We assume that the initial Lorentz factor of each shell is 
distributed uniformly in logarithmic space between log7 m j n and log7 max . The initial masses 
are assumed to be correlated with the Lorentz factors as m oc r y v . For rj — — 1 and the 
shells initially have equal energy and equal mass respectively, and for rj — 1 the shells 
initially have equal density under an assumption of the equal shell width. We assume a 
constant value I for the initial widths and the initial separations between the shells. Then, 
the initial position of the shells are Ri = 2{i — 1)/. 

The evolution of the system in time is basically same as in KPS97, but equation [| is 
used to calculate the Lorentz factor of the reflecting shells for the next time step. We follow 
the evolution of shells until there are no more collisions, i.e. until the shells are ordered 
with increasing value of the Lorentz factors. 

The conversion efficiency from the kinetic energy of the shells to radiation can be 
calculated by using the initial and final kinetic energy as (e) = 1 — Sm- 7- /Em- 7- 
where the superscripts (/) and (i) represent the initial and final values respectively. It 
depends on the model parameters {jmax/jmin, N, r),e e ,8} and on the specific realization: 
the set of random Lorentz factors assigned to each shell. For each choice of the parameters 
of the model, we have evaluated the efficiency for 100 realizations. The mean efficiency and 
its standard deviation are listed in Table 1 for the conserved equal mass case (77 = 5 = 0). 



Table. 1 
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N 


imin 


^fmax 


efficiency [%] (e e = 0.1) 


efficiency [%] (e e = 0.5) 


30 


10 2 


10 3 


9.2 ±2.3 


16.2 ±3.0 


30 


10 


10 4 


40.0 ±9.2 


67.5 ±9.3 


10 2 


10 2 


10 3 


15.1 ±1.5 


17.7 ±1.5 


10 2 


10 


10 4 


62.9 ±4.8 


72.4 ±4.3 



The efficiency approaches an asymptotic value as the Lorentz factor ratio •y m ax/lmin 
and the number of the shells iV increase. The asymptotic value depends on rj, e e and 5. The 
efficiency is plotted in fig |3| as a function of rj for 5 = 0, 0.4 or 0.8. In the range — 1 < r\ < 1, 
it is not very sensitive to 5 and peaks around r\ ~ 0, i.e. the most efficient case is the equal 
mass case. 

The efficiency is plotted as a function of e e for the conserved equal mass case (rj = 5 = 0) 
in Figure [|. It is interesting that the efficiency can be larger than e e since the energy 
released in GRB can be much larger than that in the afterglow. For instance, consider 
the case of N = 100, 7 mi „ = 10, 7 max = 10 4 , rj = 5 = and e e = 0.1, the internal shocks 
can convert about sixty percent of the kinetic energy to the radiation. The reminder is 
converted to the thermal energy by external shock. However, only the fraction e e = 0.1 of 
that is emitted as the afterglow Q, the efficiency by the external shock is only four percent. 
The GRB to the afterglow ratio for these parameters is about 15!. 

Figure ||a shows the resulting temporal structure in the equal mass case. It is a 
superposition of pulses from the elementary two shell collisions. Though numerous collisions 
take place during the evolution, the number of peaks in the profile is of order of N. Since 
all peaks have widths of the same order of magnitude, the different amplitude of the peaks 
originates mainly from the difference in the internal energy produced by the collisions. In 
Figure |||b, we plot the initial Lorentz factor as a function of time when shells were emitted 
by the source. We evaluate E int for all pairs of the shells if inner shell is faster. We assign 
E int to the ejection time of the inner shell (see Figure f)]c). It resembles well the temporal 
profile in [5|a. Therefore, despite the complicated nature of the two shells collisions, the 
observed burst closely follows the inner engine temporal profile. 

Neighboring shells collide on a time scale of 2 / y 2 l/c. The matter is moving toward 
the observer, the resulting observer time scale is l/c. On the other hand, the difference 
in observer time due to the location of the given shell within the wind is of order Nl/c. 
Then, we observe the pulses arising from the collisions according to their positions inside 
the wind. In Figure |^d we plot the shell's index against the time when a radiation from the 



: In fact, Sari 1997 have shown that also during the afterglow the fraction of energy that can be emitted 
may exceed e e . However, this will be spread over many decades of time. 
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shell is observed. Although the light curve in figure ||a is the superposition of all pulses, in 
figure |5|d we plot only the pulses which are higher than a tenth of the highest. We can see a 
clear correlation between the index and the time, the temporal profile reflects the activity 
of the source. 

The deviation from the correlation, e.g. some groups of circles in a line from left top 
to right bottom is due to that a rapid inner shell collides into the outer neighbor, and the 
boosted shell further collides into the outer. Since the masses of the shells are equal, the 
Lorentz factors of the inner and the outer shells are just switched at each collision if the 
radiation is negligible. Then, if an initial Lorentz factor of a shell is peculiarly high, the 
index of the shell which has the high Lorentz factor propagates outwards. In this sequence 
with the radiation loss, the pulse from the collisions damps quickly, then, the overall 
correlation is not destroyed. 

The temporal structures for equal energy cases (rj = —1) are plotted in figure |6|. The 
merger at each collision is assumed to split to the original masses (figure ||a) or to the 
modified masses 5 = 0.4 (figure |]c). The initial distributions of 7 are the same with that of 
figure H For 77 = —1, E int takes almost the same value for most of the collisions (KPS97), 
so most of the peaks have comparable amplitudes. The profiles still reflect the activity of 
the source well. 



6. Conclusions 

We have shown that the conversion efficiency from kinetic energy of relativistic shells 
to radiation can be close to hundred percent if the source produces shells with comparable 
masses with different Lorentz factor, especially when the logarithms of the Lorentz factors 
are distributed uniformly (see also Beloborodov 2000). It had been assume in our previous 
work that the Lorentz factors themselves were distributed uniformly, and the efficiency was 
less than forty percent. With this distribution, the most efficient case is that the source 
produces shells with comparable energy, instead of comparable masses (KPS97). 

However, this high efficiency is achieved assuming that all the internal energy is emitted 
at each collision, this is not reasonable since only electrons radiative effectively, and these 
do not have all the internal energy. After the electrons radiate, large amount of the internal 
energy remains in protons. Using a hydrodynamic simulation, it has been shown that the 
hot merger produced by a collision spreads to transform the remaining internal energy back 
to the kinetic one. A simplified description of this process is to assume that the shells 
reflect each other with a smaller relative velocity, after the collision. The difference of the 
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kinetic energy gives the radiated internal energy. 

Since the reflecting shell collide into the outer neighbor shell, the index of the shell 
which has a high Lorentz factor propagates outerwards until the high value decays by 
the radiation loss. The shell itself which has initially a high Lorentz factor might not go 
through many collisions, but its high "kinetic energy" does. Therefore, the internal shock 
process is very efficient even if the fraction of internal energy emitted at each collision is 
small. Previously, the efficiency in the case of e e < 1 had been estimated as smaller by a 
factor of e e than that in the corresponding fully radiative case. Our ultra efficient internal 
shocks scenario shows this to be a significant underestimate. Though the efficiency that we 
have estimated is bolometric efficiency, the efficiency from the kinetic energy of the shells to 
gamma-ray radiation is also high if the faction of the energy radiated in the BATSE band 
is not very small. 

Numerous collisions happened in our ultra-efficient internal shocks model, it made the 
peak width wider than in the previous internal shocks model. However, the number of main 
peaks is still almost the same as the number of shells that the source emitted. There is a 
strong correlation between the time at which we observe a pulse and the emission time of 
the corresponding shell from the source. This correlation persists even for a small e e case 
where larger number of collision happen. The temporal structure reproduces the activity of 
the source. 

We have shown that the efficiency of the internal shock process is not limited by e e , 
while that of the external shock is so. If a fraction ( of kinetic energy of an explosion is 
converted to the radiation by the internal shocks, all the remaining one is converted to the 
thermal by external shock in the afterglow stage, and a fraction e e of the thermal is emitted. 
Then, we can roughly estimate the ratio between the energy released in afterglow and that 
in GRB as ~ e e (l — C)/C Assuming e e = 0.1, the ratio is 1/10 for ( = 0.5 and decrease as ( 
increase. If the efficiency of internal shocks is indeed very large, the luminosities of GRB 
and the afterglow are expected to be anticorrelated. 

As a general argument in internal shock model, the luminosity distribution of multi 
peak bursts would be narrower than that of bursts with only few peaks (Kumar & Piran 
2000), since the number of the peaks is almost the number of the shells N as we have 
shown. We verified that the dispersion of the efficiency is proportional to 1/ v^V in our 
numerical model. 
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Fig. 1. — Fully radiative case (e e = 1): efficiency vs A . thick solid (numerical simulation: 
average of 100 random simulations with 100 equal mass shells), thin solid (analytic estimate), 
dashed ((e) = A 2 /2) 




Fig. 2. — (a) Lorentz factor 7 vs distance x from the contact surface at t — 5000 (thin) and 
6000 (thick). The Lorentz factor 7 and position x are measured in the rapid shell frame, 
and the density p is the the fluid local frame, (b) mass fractions: 7 > 5 (solid) and 7 < 5 
(dashed) 
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Fig. 5. — Temporal structure, (a) Numerical temporal structure, (b) Lorentz factor 
of a shell at the ejection vs ejection time (c) E int between neighbor boring shells vs 
ejection time (d) Index of a shell vs observed time of the radiation produced in that shell. 
N = 100, 7min = 10 2 , lmax = 10 3 ,7] = 5 = and e e = 0.1. 
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Fig. 6. — Temporal structure and correlation between the index of shell and the radiation 
observed time. The merger at each collision is assumed to split to the original masses (a 
and b) or to the modified masses with 5 = 0.4 (c and d). N — 100, r\ = —1, 7 min = 10 2 and 
Imax = 10 3 and e e = 0.1. 



